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Abstract: 

We obtain relaxation times for field theories with Lifshitz scaling and with holographic 
duals Einstein-Maxwell-Dilaton gravity theories. This is done by computing quasinormal 
modes of a bulk scalar field in the presence of Lifshitz black branes. We determine the 
relation between relaxation time and dynamical exponent z , for various values of boundary 
dimension d and operator scaling dimension. It is found that for d > z+l , at zero momenta, 
the modes are non-over damped, whereas for d < z + 1 the system is always overdamped. 
For d = z + 1 and zero momenta, we present analytical results. 
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1 Introduction 

Much effort has been invested into the understanding of the application of holography since 
its seminal papers [1-3]. Motivated by the fact that real world systems often exhibit non- 
relativistic scale invariance at critical points, rather than relativistic invariance, it is of 
interest to extend the holographic dictionary into non-relativistic context [4-8]. As a model 
of non-relativistic scaling at a critical point on the boundary we require invariance under 

t — X z t , x —^ Ax, (1.1) 

where t denotes time, x denotes spatial coordinates on the holographic boundary and z is 
called the dynamical exponent. For z = 1 the theory is Lorentz invariant. When z > 1 
the system obtains anisotropic scaling between space and time, often called Lifshitz scaling, 
which violates Lorentz boost invariance. 

The main goal of this paper is to determine the relaxation time of the d-dimensional 
boundary theory with dynamical exponent 2 at some temperature T. Relaxation occurs 
after having perturbed the system by an operator at the boundary. We consider operators 
with scaling dimension A and a spin-zero held as holographic dual in the bulk. We ob¬ 
tain the relaxation time by computing quasinormal modes in the bulk and investigate the 
dependence on d, z and A. 

To obtain Lifshitz scaling and temperature on the boundary, we consider a black brane 
solution of an Einstein-Maxwell-Dilaton (EMD) action [8, 9]. A black brane in the bulk 
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drives field excitation into dissipation. We then solve the (complexified) Klein-Gordon equa¬ 
tion of the spin-zero field in the probe limit. The resulting complex eigenvalues are called 
the quasinormal frequencies corresponding to quasinormal modes. The smallest imaginary 
part of these eigenvalues is inversely proportional to the relaxation time r of the boundary 
system. 

The quasinormal modes for z = 1 with d = 3,4,6 were studied numerically for the 
first time in [10]. An analytic solution for d = 2 , z = 1 was obtained in [11], In this current 
paper we find a generalization of this analytic solution for d = z+l for vanishing momenta. 
It reads 

= -UwT (in + A) , =* (1.2) 

where T denotes temperature and n = 0 gives the lowest lying quasinormal mode. For 
d = 3, z = 2 this equation coincides with the one found in [12]. For reviews on this topic 
we recommend [13, 14]. 

A summary of previously obtained analytic and numerical solutions of quasinormal 
modes of spin-zero fields, making usage of various bulk actions, entails: d = 2, z = 3 in 
New Massive Gravity (NMG) [15], d > 4 , z = 2 in a R 2 gravity setting [16] and d > 2 , z = 2 
in a R 3 gravity setting [17]. The case for d = 3, z = 2 has been studied in the Einstein- 
Proca-Scalar (EPS) background [6, 18], in the EMD setup [12] and in a topological black 
hole in a Einstein-Maxwell-Proca (EMP) background [19-21], All these quasinormal mode 
solutions were found to be purely imaginary. This signals that the corresponding system 
is overdamped. In the papers [12, 16] it was therefore conjectured that for (most) Lifshitz 
black holes the quasinormal modes are purely imaginary. Our conclusions will be different. 
In our numerical analysis we find that for d > z + 1 the quasinormal modes have a real 
component. However, for the case of d < z + 1 one continues to find overdamped solutions. 

A brief outline of the paper is the following. In the next section we introduce notations 
and derive the Schrodinger-like equation for quasinormal modes of a spin-zero held. In 
Section 3 we obtain the new analytic solution for quasinormal modes and analyze the 
remaining cases using numerics. Here we chart the (non-)overdamped region. Next, in 
Section 4 we present the relaxation times and their behavior versus z, d and A. Finally, 
in Section 5 we present an outlook for possible future work. 


2 Schrodinger problems for Lifshitz geometries 

The objective is to obtain relaxation times for a boundary held theory probed by an oper¬ 
ator dual to a spin-zero held in the bulk. In order to pursue this goal we need to compute 
quasinormal modes of a spin-zero held in an asymptotically Lifshitz black brane back¬ 
ground. This requires solving the equation of motion of a massive scalar held in the EMD 
background, in the probe limit and subjected to appropriate boundary conditions. 
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2.1 Lifshitz brane 


The metric line element of a black brane exhibiting Lifshitz scaling (1.1) on the boundary 
can be expressed as 


ds 2 = 


1 


7"2^y2 


dr 2 — V 2 {r)r 2z dt 2 + r 2 dx 2 l _ 1 , 


( 2 . 1 ) 


where r is the radial bulk coordinate, which under Lifshitz symmetry scales as r —> A r. 
The limit r —> oo corresponds to the boundary. The Lifshitz radius, a generalization of 
the Anti-de Sitter radius, is put to unity in this paper. The blackening factor in the EMD 
setup 1 is [8] 

f r h \d+z-i 


V z = 1- ( — 
r 


) 


where denotes the horizon. The temperature of the black brane is given by 

47tT = (d + z — 1 )rf l . 


(2.3) 


(2.4) 


To find the quasinormal modes we need to compute the tortoise coordinate r* , which 
characterizes the radial null curves obeying t = ±r* + constant. We consequently demand 


dr * = r z + ly2 dr - 


The general solution, for d + z > 1, is 


r* = -2 Fi 

z 


1, 


d + z — 1 


; i + 


d+z—l\rJ 


r _ h \ d + z - 1 


(2.5) 


( 2 . 6 ) 


in terms of the hypergeometric function of the second kind. Without any loss of generality 
we fix the integration constant to be zero. For the special case of d = z + 1 it becomes 


1 


*\d=z+l 


2 zr z h 


log 


1 + ^ 


1 — 


Next we define Eddington-Finkelstein coordinates v and u as 

v = t + r* , u = t — r* , 

and we require infalling boundary conditions, for a held (f> near the horizon to be 

<J>(r -> r h ) ~ e~ iujv = e~ iut - iujr *. 

This condition causes dissipation and complexihes the held 0. 


(2.7) 


( 2 . 8 ) 


(2.9) 


1 For different bulk fields, such as e.g. the EPS setup in d = 3 and 2 = 2, one has [6, 18] 


T e 2 ps = 1 - . 
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( 2 . 2 ) 













2.2 Quasinormal modes of a scalar probe 

We consider the equation of motion of a massive scalar field 0 on the background of (2.1) 

□i f> = m 2 (j). (2.10) 


Assuming the probe limit we ignore back reaction on the metric. We require boundary 
conditions such that at the horizon the held i j) has to be falling into the black brane and 
that at radial infinity the held goes to zero, <j)(r —> oo) = 0. This system has a dissipative 
nature since over time more and more held excitations will fall into the black brane while 
the excitations at radial infinity are rehected. There is no incoming hux from the horizon 
nor from the boundary. Quasinormal modes are the discrete values of energy u for which 
the equation of motion of <j) is satished, taking into account these boundary conditions. 
Making a plane wave Ansatz and a radial rescaling 

<j> = ^{r)e-^ t+its , (2.11) 

where k denotes momentum, the Klein-Gordon equation in (2.10) can be written as a 
one-dimensional time independent Schrodinger-like equation 

[d 2 + u? - V(r)] 0 = 0, d*=r z+1 V 2 d r . (2.12) 


The corresponding potential is 


Y(r) = r 2z V 2 (r) ( -^ + m 2 + 


k 2 


(d — l)(d + 2z — 1) (d— 1) 


+ 


-) 


d-\-z —1 


(2.13) 


Asymptotically, near the boundary, the solution to the Schrodinger-like equation is given 

by 

cj)(r) ~ Ar~ A ~ + Br~ A +, A± = — + ~ -- ± \ j + ^ ^ + m 2 , (2-14) 

where A and B are independent of r. Applying the regular recipe for holography [2, 3], we 
consider the term containing A = A + to be the normalizable mode and hence the scaling 
dimension of the operator 0\ dual to 4> is A. Considering the Klein-Gordon inner product 
we find that for this mode to be normalizable we have to constrain the scalar mass [22] 

/ drr d ~ z ~ 2 \Br~ A \ 2 < oo , m 2 > — ^ ^ =m 2 BF , (2-15) 

where the contribution is due to a volume factor resulting from taking a space-like 

slice as integration area. The m 2 BF is the minimal value above which m 2 has to remain. 
This is the Breitenlohner-Freedman (BF) bound. We require 


A(ui, k , A, d,z,T) = 0 , 


(2.16) 


in order to have no incoming or outgoing flux at radial infinity. 








The BF bound (2.15) translates into the following bound on the scaling dimension 


A = 


d + z — 1 
2 


+ 



2 

+ m 2 > 


d + z — 1 
2 


(2.17) 


Marginal operators in a theory with Lifshitz scaling have A marg i na i = d + z — 1 , so we 
summarize 


d + z — 1 


< A 


relevant 


< A 


marginal 


= d + z — 1<A 


irrelevant 


(2.18) 


which can be visualized as in Figure 1. A quasinormal mode is defined by solving the 



Figure 1: Visualization of (2.18). Exactly at the d = 1 — z + A line the operator is 
marginal. On and above the line d = 1 — z + 2A the BF bound is violated. 


Schrodinger-like equation (2.12) supplemented with boundary condition (2.9) at the horizon 
and boundary condition (2.16) at radial infinity. Quasinormal modes u n are then obtained 
from 


A(u, k, A, d,z,T) = 0, 


ui n = ^Re{n,k,A,d,z,T) - iu}i m (n, k, A, d, z, T), (2.19) 


where integer n labels the overtone number of the quasinormal mode. By definition the 
n = 0 gives the imaginary component closest to zero. 

When oj\ m is positive (in our conventions) it implies stability of the gravitational back¬ 
ground under scalar perturbations. We can verify that uj\ m is always positive due to the 
fact that the potential is real, positive and strictly increasing. This conclusion follows from 
a reasoning similar to the one given in [10], but adapted to Lifshitz scaling. One starts by 
introducing cj) = e~ lur *(j)(r) to (2.12) in order to obtain 


r z+1 V 2 { r )d 2 r + [d r (r z+1 V 2 (r)) - 2 iu] d r 


V(r) 

r z+ly2 


= 0 . 


( 2 . 20 ) 


Multiplying this equation by (f>* and integrating over r from 77 , to 00 , after integration by 
parts, yields 



r z+1 V 2 (r)\d r 4>\ 2 + 2ioj(f)*d r (f) + 


V(r) 

r z+iv 2 (r) 



= 0 . 


( 2 . 21 ) 















Taking the imaginary part of (2.21), applying integration by parts and inserting the result 
back into (2.21) results in 


/*oo 

/ dr 
Jrh 


r z+ 1 V 2 (r)\d r (f >\ 2 + 


V(r) 


r z+iV 2 (r) 


I 


\u\ 2 \4>( r h)\" 

Imu 


( 2 . 22 ) 


where the left-hand side is ensured to remain positive with the potential under consideration 
in this paper. This guarantees a negative imaginary value of to. 

The connection between relaxation time and quasinormal modes follows from 




(2.23) 


which violates unitarity due to dissipation. The relaxation time t is defined as 


r = 1/wim . (2.24) 

Notice that the real part of the quasinormal mode is referred to as the dispersion relation. 
For purely imaginary quasinormal modes we thus speak of the corresponding system being 
overdamped. 


3 Obtaining quasinormal modes 

Finding the quasinormal modes by solving the Schrodinger-like equation (2.12) supple¬ 
mented with boundary condition (2.9) at the horizon boundary condition (2.19) at radial 
infinity is mostly done numerically in this paper. The case d = 2 + 1, however, can be 
treated analytically for vanishing momenta, so we can consider this case separately. 

3.1 Analytic solutions, d = z + 1 

For the EMD background, an analytical solution of quasinormal modes is known d = 3, 
2 = 2 with k = 0 [12]. We now determine additional analytical solutions for the general 
case of d = 2 + 1 with k = 0. 

Why this case is special from an analytic point of view can be understood from observing 
two simplifications which occur to (2.12) in this case. Firstly, the potential only depends 
on terms containing r 2z . This simplifies the r-dependence of the potential V tremendously. 
Secondly, r* depends only on r z and can be analytically inverted to r 2 (r*), which enables 
one to directly express the potential V in terms of r* . In other cases one is unable to 
accomplish this or it simply yields a much more complex expression. 

However, in order to compute this analytic expression it is convenient to switch back 
from 4> to (j) , require d = z + 1, and write (2.12) in the form 

r 2z V 2 (r)m 2 + J 2 — r 2z V 2 4> + d 2 (f> + r z+1 V 2 {r)-d if (j) = 0 . (3.1) 

We make use of the substitution 

r h \ 2 






in order to obtain 


UJ 


l y z — r h Z 2 + k 2 y] (1 — y z ) (z — 1) + y z , „ , , 

y h 1 h yj v y ' ^ - 4 ^ , ’ y + 42/0" = 0, (3.3) 


r 2z(i _ yZ y y 

which, when assuming k = 0, is solved by 
(, b=Ciy ~ 2 ~ (1 - y z )P 2 F\ 

+ c 2 y~^ (l - y z Y 2 F 1 


l ~y z 






(3.4) 


In the last expression c\, c 2 are some constants and 

0= ” 


AnT 


(3.5) 


The deltas are the same as given in (2.14). We continue by requiring both boundary 
conditions. We start with the boundary condition at r —> 00 , which corresponds to y —> 0. 
We want Dirichlet boundary conditions to hold here. Taking the limit one obtains 


4>(y —> 0) ~ ciy 2 + c 2 y 2 ' ~ c\r A - + c 2 r ^ + , 


„-a 4 


(3.6) 


thus we put ci = 0 in order to match the required boundary behavior from (2.16). When 
taking the limit we made use of the identity 


2 Fi[a, b ; c; 0] = 1 for all a, b , and c non-zero, 


(3.7) 


which can be found in e.g. [23]. Next we apply the boundary condition at the horizon 
r = rh , which corresponds to y = 1. We can relate the behavior at y = 1 to the behavior 
at y = 0 using 


(3.8) 


_ r(c)T(c — 2 a) 

2 F\[a, a; c; z = ——-— 2^1 a, a; 2a - c+ 1; 1 - z 

I (c — ay 

+ (1 _ z)C -2aI>)r(2a-c) 2i r 1 [ c _ 0jC _ a;c _ 2a + 1;1 _ (Z ] j 
I (a) 

where T denotes the gamma function. Using (3.7) we arrive at 

+ = iim (1 - y-f / r(A + / z) r(-2ffl _ „ r(A + /s)r(2ft i 

y^ lr(A + /(2 z)-Df r(A + /(22) + /j) 2 J 

Now we have to identify the ingoing modes as written down in (2.9). From (2.7) we collect 

1 


(3.9) 


*\d=z+l 47tT 

such that we are able to identify 


[log (1 — y z ) — 2 log (1 + y z )\ , 


(3.10) 


e~ iut (1 - y 2 ) ±/3 = e- iut e ±iuj l °^T ] = + , 
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(3.11) 















which leaves us to the conclusion that we have to put the term with the positive power of 
/3 to zero. This is accomplished by requiring 


r(A+/(2«)- j 8) °’ ^ 2irT 



(3.12) 


where n is a positive integer. For z = 2 the result restores the finding in [12]. For z = 1 
we obtain the case of [11]. This solution holds for any d > 1 and at all times exhibits 
overdamped behavior. It is concluded from (2.24) that 

Td=z+ 1 = 2vrTA ’ ( 3 ' 13 ^ 

which implies that higher anisotropy corresponds to a longer relaxation time. From the 
numerics it is observed that this holds as long as z > d — 1. 


3.2 Numerical solutions, d / z + 1 

Initially we started by solving the Schrodinger-like equation (2.12) using the method de¬ 
scribed in [10]. This method, in essence, solves the Schrodinger-like equation by using a 
power series Ansatz. Applying the boundary conditions, results in a recursive relation be¬ 
tween different terms in the power series Ansatz. This recursive relation is too cumbersome 
to be handled analytically. We thus have to resort to numerical methods. 

The time it takes for the power series to converge numerically, is found to increase 
dramatically for z > 1. It was already noted by the authors of [10] themselves that 
increasing d increases computation time. This is connected to the fact that increasing d, 
increases the power of r in the potential, yielding a more complex recursive relation, which 
amplifies the time it takes to converge. Noting that when increasing z that the power of r 
in the potential gets larger too, explains the rise of computation time. 

To decrease computational time we adopted the Improved Asymptotic Iteration Method, 
as described in [24], This method, however making use of a recursive structure as well, relies 
on (as opposed to the previously mentioned algorithm) the observation that 

x"( x ) = A 0 (x)x'(x) + s 0 (x)x(x ), => X (n+2) ( x ) = K( x )x'( x ) + s n (x)x(x), (3.14) 


where the superscripted (n + 2) denotes the order of derivation. The Ao and so are poly¬ 
nomials to which A n and s n are related in a recursive fashion including various orders of 
derivatives as well. From the ratio of the (n + 3) th and (n + 2) th derivatives, one shows 


^log (x (n+2) ( 


x = 


An+i (x'(z) + ££^x(®)) 


A» 


\X'( X ) 


+ 


Sn{x) 
A n ( x ) 


X(x 


(3.15) 


Now we introduce the asymptotic aspect of the method. If for some sufficiently large n 

^ ^ 71+1 


Sn 

X r , 


^n+1 


(3.16) 


one can solve (3.15) and, by plugging this back into (3.14), find a solution for x( x ) ■ This 
approach is called the Asymptotic Iteration Method and was originally developed by [25]. 


- 8 - 









The Improved Asymptotic Iteration Method entails, as modification to the original ap¬ 
proach, some convenient power series expansions of s n and A n in order to simplify their 
respective recursive relations. 

To put this algorithm to our use we choose to rescale the coordinate r into a dimen¬ 
sionless parameter with a finite range in the following way 

x = l~ — . (3.17) 

r 

We rewrite equation (2.12) into 


I w „ 2 ) «*) - 0, 


(1 — x) 2 h(x ) 


=Q(a;) 


'[(l-x) 2 /i(x)]^ 

Eii(l) 


(3.18) 


where /i(r) = r z+1 V 2 (r) and the accents denote derivatives with respect to x. We aim to 
scale out the behavior near the horizon and the boundary. This is done by employing the 
scaling 


(j>(x) = x A (l — x) B x{ x ) i A = 


— ILO 


(d + z — l)r); 


B =2 + 


m 2 + 


(d + z- l) 2 


(3.19) 


Using that 

_=nu_ 

d x \x A (l — x ) B ] = x A (l — x) B \_Ax ~ 1 — B( 1 — x) _1 ] , (3.20) 

we rewrite (3.18) as 

x"(x) + [2 F(x) + Q{x)\ x'(x) + [F 2 (x) + d x F(x ) + <5(x)i ? (.x) + i?(x)] x( x ) = 0 . (3.21) 

s ---" '-V-" 

=Ao =so 


Now that the form of (3.14) is obtained, we can in principle construct A n and s n for any 
n . The quasinormal modes oj are obtained by requiring and solving 


s n (x,u)X n+1 (x,uj) - s n+1 (x,u)\ n (x,uj) « 0 , 


(3.22) 


for some large value of n. Typically we use ~ 30 to obtain certainty up to at least two 
decimals. In our case for z > 1 we obtained quicker convergence with this algorithm, than 
when using the first mentioned algorithm. We reproduced the results from [10] for the case 
of a black brane in Table 1 in the appendix as a check. 


3.3 Analysis of numerical output 

We focus on d = 2, 3,4 because of their relevance in real world systems. We choose momen¬ 
tum to be vanishing and leave the k ^ 0 to future work. For clearness we restrict ourselves 
to A = 3,4, 5.5 and values of z for 1 through 6. To get some intuition for the behavior 
of the quasinormal modes in the complex plane, when varying z , we present a cartoon in 
Figure 2. 
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Figure 2: The hollow dots denote the location of the quasinormal modes when z = 1 . 
The n denotes the overtone number. When increasing z, the location of quasinormal mode 
will follow the arc towards the vertical axis. At z = d — 1 the mode, for the first time, 
hits the vertical axis. From there on, for any z > d — 1, the quasinormal mode remains 
somewhere on the vertical axis. We stress that remaining on the vertical axis is because of 
the vanishing real part when z > d — 1 and corresponds to overdamped systems. 



Figure 3: The real part of the computed quasinormal modes. For d = 2 the value is zero 
everywhere, as can be read off from Table 2 in the appendix. Diamonds correspond to 
A = 5.5, triangles correspond to A = 4 and the dots correspond to A = 3. 


In Figure 3 we plot the real part of the quasinormal modes versus z. A sample of 
the data can also be found in Table 2 in the appendix. In the region d < z + 1 we find 
that the system is overdamped. For d > z + 1 we find a non-zero real part. This result 
is interesting when taking into account the conjectures of [12, 16] , which state that for 
(most) Lifshitz black holes the quasinormal modes are purely imaginary. Moreover, finding 
non-overdamped cases for d > z + 1 contrasts the overdamped cases found for d > 4, 
z = 2 in a R 2 gravity setting [16] and d > 2, z = 2 in a R 3 gravity setting [17]. We 
present a cartoon of the qualitative structure of the overtones of the quasinormal modes 
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in the complex plane in Figure 4. We summarize Figures 1 and 4 in Figure 5 for clarity. 



Figure 4: The line given by d = z + 1 divides quasinormal modes in being overdamped or 
non-overdamped. Both plots within the separated regions show the qualitative behavior of 
the overtones of the quasinormal modes in those regions. Exactly at and underneath the 
line d = z + 1 the system is overdamped. 

Notice that for fixed A , d and increasing z , the real part is strictly decreasing, until it hits 



Figure 5: The information from Figures 1 and 4 is summarized. Exactly at the d = 1—£+A 
line the operator is marginal. On and above the line d = 1 — z+2A the BF bound is violated. 
Exactly at and underneath the line d = z + 1 the system is overdamped. 

d = z + 1. When fixing z we have the relation that higher A corresponds to a greater real 
part in the non-overdamped region. 

4 Relaxation times 

In Figure 6 we present the relaxation time versus z . For the relaxation times the point 
d = z + 1 leaves its footprint as well as in the real part, by separating different behaviors on 
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Figure 6: Results of the relaxation times r, numerically computed from quasinormal 
modes. Diamonds correspond to A = 5.5, triangles correspond to A = 4 and the dots 
correspond to A = 3. A star denotes that the operator at that point is marginal. The 
curve stops on the right when the BF bound is violated. The highlighted points denote 
when the operator is relevant, otherwise the point is irrelevant. However, around the region 
d = z + 1 we did not put any highlights for readability’s sake. 
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either side of this point. A sample of the data points is given in Table 2 in the appendix. 
Features which we gather from Figure 6 are: 

• Higher scaling dimension A of an operator corresponds to a lower relaxation time. 

• Increasing z , the amount of anisotropy, corresponds to a higher relaxation time r 
when d < z + 1 . 

• Overdampedness is independent of the dimension of the operator. 

5 Outlook 

The results in this paper were obtained from the holographic point of view. It would be 
interesting if, in the spirit of [11], one could reproduce these results directly from the field 
theory side. In particular it would be interesting to understand the difference between the 
regimes d < z + 1 and d > z + 1. In context of this we make the following observation. 
Consider a free d-dimensional theory with Lifshitz scaling and a dispersion 

u~\k\ z . (5.1) 

The number of states f l(k) up to momentum \k\ is obtained by taking a spherical fc-space 
volume and dividing it by the volume occupied per allowed state, resulting in 

n(k) ~ k d ~ l . (5.2) 

The density of states D(u>) is defined as 

dQ,(k(u)) dfl(k) dk(uj ) ^2 1 -* <M*+i) 

D(ui) = ---= —-- - — ~ u; * u; * = co * . (5-3) 

doj dk duj 

Notice that there is a qualitatively different behavior between the regimes d < z + 1 and 
d > z + 1. For d < z + 1, the density of states decreases with energy, whereas for d > z + 1 
it increases with energy. Perhaps this behavior is related to the non-overdamped and 
overdamped phases after adding interactions. We leave this for future research. 
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Appendix 

In Table 1 we present data in order to compare to the results of [10]. Table 2 contains a 
sample of the results used in this paper. 



d = 

= 3 

d = 

= 4 

d = 

= 6 

Th 

MIm 

URe 

MIm 

URe 

WIm 

URe 

100 

266.38 

184.94 

274.66 

311.94 

261.24 

500.74 

50 

133.19 

92.47 

137.33 

155.97 

130.62 

250.37 

10 

26.63 

18.49 

27.46 

31.19 

26.12 

50.07 

5 

13.31 

9.24 

13.73 

15.59 

13.06 

25.03 

1 

2.66 

1.84 

2.74 

3.11 

2.61 

5.00 


Table 1: Results which can be compared to [10]. Notice that this corresponds to putting 
the mass m to zero and z = 1. 
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A 

= 3 





A = 

5.5 




d = 

= 2 

d~- 

= 3 

d = 

= 4 

d = 

= 2 

d = 

3 

d = 

= 4 

z 

Re (wo) 

AtyT 

47tTt 

Re (wo) 
47 tT 

47tTt 

Re (wo) 
4 ttT 

47tTt 

Re (wo) 
AtyT 

47tTt 

Re(wo) 

AnT 

47tTt 

Re (wo) 
4 ttT 

47 tTt 

1 

0 

2/3 

0.61 

1.12 

0.54 

2.27 

0 

4/11 

1.14 

0.53 

1.13 

0.93 

4/3 

0.00 

1.03 

0.39 

1.15 

0.45 

2.01 

0.00 

0.54 

0.66 

0.57 

0.87 

0.87 

5/3 

0.00 

1.30 

0.21 

1.23 

0.36 

1.91 

0.00 

0.69 

0.32 

0.63 

0.64 

0.87 

2 

0.00 

1.55 

0 

4/3 

0.28 

1.89 

0.00 

0.82 

0 

8/11 

0.46 

0.90 

7/3 

0.00 

1.80 

0.00 

1.83 

0.204 

1.90 

0.00 

0.95 

0.00 

0.96 

0.31 

0.95 

8/3 

0.00 

2.04 

0.00 

2.15 

0.13 

1.94 

0.00 

1.09 

0.00 

1.13 

0.18 

1.02 

3 

0.00 

2.28 

0.00 

2.44 



0.00 

1.21 

0.00 

1.28 

0 

12/11 

10/3 

0.00 

2.52 

0.00 

2.72 



0.00 

1.34 

0.00 

1.42 

0.00 

1.38 

11/3 

0.00 

2.76 

0.00 

2.99 



0.00 

1.47 

0.00 

1.56 

0.00 

1.57 

4 

0.00 

2.99 





0.00 

1.60 

0.00 

1.70 

0.00 

1.73 

13/3 

0.00 

3.23 





0.00 

1.72 

0.00 

1.84 

0.00 

1.89 

14/3 

0.00 

3.46 





0.00 

1.85 

0.00 

1.97 

0.00 

2.04 

5 







0.00 

1.98 

0.00 

2.11 

0.00 

2.19 

16/3 







0.00 

2.10 

0.00 

2.24 

0.00 

2.33 

17/3 







0.00 

2.23 

0.00 

2.37 

0.00 

2.48 

6 







0.00 

2.35 

0.00 

2.50 

0.00 

2.63 


Table 2: Parts of the results which are used to plot Figure 3 and 6. Empty rows signal 
violation of the BF bound. 
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